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ABSTRACT 

We investigate several fundamental properties of z ~ 4 Lyman-break galaxies 
by comparing observations with the predictions of a semi-analytic model based 
on the Cold Dark Matter theory of hierarchical structure formation. We use a 
sample of -B435-dropouts from the Great Observatories Origins Deep Survey, and 
complement the ACS optical -B435, Veoe, "^775, and 2:850 data with the VLT ISAAC 
J, if, and Kg observations. We extract -B435-dropouts from our semi-analytic 
mock catalog using the same color criteria and magnitude limits that were applied 
to the observed sample. We find that the ^775 — Kg colors of the model-derived 
and observed 5435-dropouts are in good agreement. However, we find that the 
^775 — ^850 colors differ significantly, indicating perhaps that either too little dust 
or an incorrect extinction curve have been used. Motivated by the reasonably 
good agreement between the model and observed data we present predictions for 
the stellar masses, star formation rates, and ages for the 2; ~ 4 Lyman-break 
sample. We find that according to our model the color selection criteria used 
to select our 2; ~ 4 sample surveys 67% of all galaxies at this epoch down to 
2:850 < 26.5. We find that our model predicts a ~ 40% mass build-up between 
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the z ~ 4 and 2; ~ 3 epochs for the UV rest-frame L* galaxies. Furthermore, 
according to our model, at least 50% of the total stellar mass resides in relatively 
massive UV-faint objects that fall below our observational detection limit. 

Subject headings: galaxies: evolution — galaxies: fundamental parameters — 
galaxies: high-redshift — galaxies: photometry — cosmology: observations — 
cosmology: theory 

1. Introduction 

The " Lyman-brcak" color selection technique has been shown to be a highly effective 
means of selecting galaxies at high redshift z >2 (Steidel et al 1999; Madau et al. 1996). 
Fairly large samples of Lyman- break galaxies (LBGs) in the redshift range 2 < 2; < 3.5 have 
been compiled, and the characteristics of these objects have been studied with much vigor 
over the past several years. Optical and NIR photometric colors have been used to ascertain 
various attributes of these high-redshift objects. For example, colors spanning the rest- 
frame Balmer break (A ~ 4000 A) have been used to constrain the range of possible stellar 
ages and masses of these galaxies through application of maximum-likelihood techniques to 
synthesized stellar population models (Papovich, Dickinson, & Ferguson 2001; Shapley et al. 
2001). 

As well, many researchers have compared theoretical models of galaxy evolution with 
various aspects of the observations (Wcchslcr et al. 2001; Baugh et al. 1998). Somcrvillc, 
Primack, & Faber (2001, hereafter SPF2001) have used semi-analytic models to perform 
detailed comparative studies on the 2; ~ 3 sample. In SPF2001 it was shown that the 
luminosity function and rest-frame UV and UV-optical colors of 2; ~ 3 LBGs could be 
reproduced well by the models, implying that the ages and dust contents of the model 
galaxies were consistent with the observations. 

In this letter, we compare the colors of observed 5435-dropouts, selected from ACS 
imaging from the Great Observatories Origins Deep Survey (GOODS) and supplemented 
with deep NIR (JHRg) data from ISAAC on the VLT, with the colors of model-derived 
-B435-dropouts. We then present model predictions for physical properties of these objects, 
such as stellar masses, ages, and star formation rates. We also investigate selection effects 
on the stellar mass content for the redshift that we sample, and we look at the stellar mass 
build-up between the 2; ~ 4 and 2; ~ 3 epochs. 

All quoted magnitudes throughout the letter are in the AB magnitude system (Oke 
1974). 
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2. Observations and Measurements 

We use three epochs of HST ACS observations of the Chandra Deep Field South (CDF- 
S) obtained as part of the GOODS program. The data consist of image mosaics in all 
four ACS bands S435, V^soe, ^775, and ^§50 (the ACS F435W, F606W, F775W, F850LP filters, 
respectively) spanning wavelength range of 0.4 to 1 iim, and covering approximately a field of 
view of 10' X 16'. The specifics of data acquisition and reduction, as well as object detection 
and photometry can be found in Giavahsco et al. (2003a). 

We supplement our AGS data set with the VLT ISAAC near-infrared imaging. The 
NIR data consists of the J, H, and Kg bands that extend our wavelength coverage out to 
A ~ 4400 A(rest-frame 2; ~ 4). We combine the AGS and ISAAC data using a method that 
optimally matches the relative photometry between these bands (Papovich 2002; Papovich 
et al. 2003). 

3. Theoretical Model 

We make use of a semi-analytic model, based on the hierarchical structure formation 
paradigm in a AGDM cosmology (Qa = 0.7, fi^ = 0.3, h = Hq/IOO kms-^Mpc"^ = 0.7). 
The model treats the formation of structure via a hierarchical 'merger tree', and includes 
a treatment of gas cooling, star formation, supernova feedback and chemical enrichment, 
galaxy mergers, stellar populations and dust (see Somerville & Primack (1999) and SPF2001 
for details). We follow halo merger histories down to a circular velocity of Vc = 30kms^^, 
as gas collapse and star formation in smaller halos is assumed to be suppressed by the 
presence of a photoionizing background (see Somerville 2002, and references therein). We 
use the multi-metaUicity stellar SEDs from the STARDUST models (Devriendt, Guiderdoni, 
& Sadat 1999), with a Kennicutt IMF. Dust extinction is modeled in a similar manner as in 
SPF2001. We use this model to produce a 'mock-GOODS' catalog with the same geometry, 
sky area, filter passbands, etc. as the real GOODS. 

4. Galaxy Sample 

The S435 -dropout sample was chosen using the color selection criteria described in de- 
tail by Giavalisco et al. (2003b). The sample was refined by applying detection limits at 
> 5 for the AGS data set, and another one at (S/N)^'^ > 1 for the matched NIR 
set. These signal-to-noise restrictions limited our sample to galaxies with 2:850 < 26.5, which 
helped curtail the presence of spurious detections. The relatively faint Ks cut-off was ne- 
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cessitated by the desire for a large enough sample for our study and did not significantly 
affect our conclusions. After further visual inspection, we chose 136 color-selected galaxies 
for analysis. 

Using the same color selection criteria that were applied to the observed sample we 
produced a model-derived 5435-dropout sample. We then applied a magnitude cut-off of 
z < 26.5 in order to comply with the observational detection limit. Before applying the color 
selection criteria we incorporated simulated observational scatter into our model-derived 
photometric catalogs. The observational scatter was drawn from a Gaussian distribution 
with the typical signal-to-noise values found in the observed CDF-S data set. 

Due to the dearth of spectroscopic data, empirical Monte Carlo simulations were per- 
formed to estimate the redshift distribution of the observed 5435-dropout sample. These sim- 
ulations were based on artificial LBGs distributed over a wide redshift range (2.5 < z < 8) 
with assumed distribution functions of UV luminosity, SED, morphology, and size, adjusted 
to match the colors of observed S435 -dropouts observed at z ~ 4. Reader is urged to see 
Giavalisco et al. (2003b) for a detailed discussion on these simulations. The redshift distri- 
bution of the simulated color-selected sample was found to have a mean value of z ~ 3.78 
with a standard deviation of ±0.34. 

One of the interesting tests that we can perform is to test for the incompleteness of our 
color-selected sample. According to the Monte Carlo simulations, 73% of all the simulated 
-B435-dropout galaxies down to 2^850 < 26.5 and in the interval 3.44 < z < 4.12 are recovered 
using our color selection criteria. When we apply the color criteria to our model, we select 
67% of all model galaxies in the same 3.44 < 2; < 4.12 redshift range, down to the same 
limiting magnitude of zg^Q < 26.5. This implies that the simulations performed by Giavalisco 
et al. (2003b) and the semi-analytic model show concordant incompleteness estimates with 
respect to the i?435-dropout selection technique in the above redshift range and down to our 
detection limit. 



5. Galaxy Properties 

The 2775 — Kg vs. 2775 — 2:850 (1550 - 4400 vs. 1550 - 1700 A, rest-frame z ^ 4) colors 
for the observed and model-derived i?435-dropout samples are presented in Figures la and 
lb, respectively. A Kolmogorov-Smirnov test reveals a 19% likelihood that the ^775 — Kg 
colors are drawn from the same underlying distribution for the two S435-dropout samples. 
A probability of 2 X 10~^ is obtained for the corresponding ^775 — 2;85o colors. This indicates 
a relatively good agreement in the i775 — Kg color distribution between the two samples, but 
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a poor correlation in the i/ys — 2^850 colors. 

The ^775 — 2:850 colors probe the slope of the UV continuum, which is believed to be 
primarily an indicator of internal dust content in young stellar populations (e.g. Meurer, 
Heckman, & Calzetti 1999). In our model we used dust normalized against the z — 
data and a Galactic extinction curve. The apparent disparity could be fixed by employing 
a different extinction curve (Calzetti 1997), or including the expected dependence of the 
extinction on the age of the stellar populaton (Chariot & Fall 2000). 

Aside from the photometric analysis, we compared the luminosity function from our 
model with the one for G-dropouts from Steidel et al (1999) and B-dropouts from HDF-N 
to ensure proper number counts. We found good agreement, which is not surprising since 
our model is very similar to the one that was used (and optimized) in SPF2001. In fact, the 
reader is urged to refer to SPF2001 for a detailed discussion on this topic. 

Several researchers (e.g. Papovich, Dickinson, & Ferguson 2001; Shapley et al. 2001) have 
used simple parameterized star formation histories and stellar population synthesis models 
to estimate physical properties of LBG's at 2; ~ 3 using the optical-NIR photometry as a 
constraint on the star formation history. Here, as we have shown that our model reproduces 
(modulo dust) the observed color distribution for objects selected via their ^435 — Vqqq and 
Vfsoe — ^775 colors and corresponding magnitude limits, we argue on similar grounds that our 
model should reproduce the statistical distribution of underlying stellar ages and masses of 
the observed S435-dropout population. Thus we can use our models to obtain estimates for 
some of these quantities. 

In Figure 2 we show the stellar masses of our color-selected model galaxies. We note that 
the stellar masses range from 10^ to lO^'^ h~^ M©, which is roughly two orders of magnitude 
less than the stellar masses of the present day spirals and ellipticals — this indicates 
that, as with the ^ ~ 3 population, the ^ ~ 4 Lyman-break galaxies are not the fully 
assembled progenitors of the present-day L > L^, galaxies (Giavalisco, Steidel, & Macchetto 
1996; Steidel et al. 1996a), and that several generations of merging events must take place 
between z ^ A and the present epoch. The median mass of Log[Mstar] ~ 9.26 (Mq h~^) is 0.5 
dex less than for the 2; ~ 3 galaxies studied in SPF2001. This implies a stellar mass build-up 
between the two epochs. To further explore the last point we have looked at the mean stellar 
masses of all L* galaxies measured in rest-frame UV and predicted by our model at 2; ~ 3 
and 2; ~ 4. We find that at 2; ~ 3, with m* = 24.358 (UV rest-frame), we get a mean 
value of Log[Mstar] = 9.74 (Mq h~^) in a m* ± 0.5 magnitude interval, and at 2; ~ 4, with 
m* = 24.998 (UV rest-frame), we get a mean value of LogfMgtar] = 9.58 {Mq h~^), again 
in a m* ± 0.5 magnitude interval. This corresponds to a mass build-up of approximately 
~ 40%. This type of mass build-up between 2; ~ 4 and 2; ~ 3 is similar to that inferred from 
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the Papovich et al. (2003) comparison of volume- averaged SEDs of observed Lyman-break 
galaxies. 

Our model also provides us with the star formation rates (SFR) and ages of our color- 
selected model galaxies. Our model results tell us that the values of these two quantities 
for the 2; ~ 4 sample are very similar to the ones found for the 2; ~ 3 sample studied by 
SPF2001. We observe a branch of actively star-forming galaxies, and a broader branch of 
galaxies with lower SFR, corresponding to bursting and fading galaxies. The star formation 
rates for the brightest galaxies in the model approach 100 Mq yr~^, similar to what has been 
found for 2; ~ 3 Lyman-break galaxies studied by SPF2001. 

We have also looked at the distribution of stellar-mass-weighted mean stellar ages. We 
found that the distribution is very broad and skewed toward ages of less than 300 Myr, with 
a median age of 200 Myr, and a peak at about 100 Myr. There is, however, a tail of objects 
with older stellar populations, reaching a mean age of ~ 1 Gyr — close to the age of the 
universe at that redshift in our adopted cosmology (~ 1.5 Gyr). 

Drawing on the stellar mass distribution results (Figure 2) we plot in Figure 3 the mass 
distribution of all (i.e. not just color-selected) galaxies from our mock catalog, limited to 
galaxies with Log[Mstar] > 9.26 (M© h~^) (the median value from Figure 2) and spanning the 
redshift range of 3.44 < z < 4.12. As was mentioned in §4, we select 67% of model galaxies 
in this redshift range with our color selection criteria, down to a limiting magnitude of 
2^850 < 26.5. This corresponds to 62% of the total stellar mass available in that redshift range, 
down to that magnitude limit. Hence, with our color criteria, we 'observe' the majority of the 
stellar mass that resides in galaxies brighter than 2:850 < 26.5 and spanning 3.44 < z < 4.12. 
Figure 2 indicates that 50% of our color-selected objects lie above Log[Mstar] > 9.26 (M0 h~^), 
we find though that only 50% of all model galaxies more massive than Log[Mstar] > 9.26 (Mq 
h^^) are brighter than 2:850 < 26.5. So for model-derived objects residing in 3.44 < z < 4.12 
we only sample 50% of the available stellar mass down to 2:359 < 26.5. The rest of the stellar 
mass resides in relatively massive UV-faint galaxies. The color selection incompleteness down 
to 2:350 < 26.5 and the mass contained in the UV-faint objects with magnitudes 2:850 > 26.5 
are two effects that conspire to severely hmit the amount of mass selected with optical surveys 
(Franx et al. 2003). Given how well our model predicts the observed colors (modulo dust), 
number counts, and other properties, this points to a substantial deficit in optically-selected 
galaxies, and by extension, the total stellar mass. Since our model catalog was limited to 
galaxies with 2:350 < 28.0 values, the 'unseen' mass fraction estimate should be taken as a 
lower limit. 
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6. Summary and Conclusions 

We performed a comparative analysis of the color-selected 5435-dropout samples taken 
from GOODS observations and a variant of a ACDM-based semi-analytic model. We found 
that the color selection technique used to obtain the -B435-dropout sample does a relatively 
good job in selecting a complete census of galaxies spanning z ~ 3.78 ± 0.34 and down to 
^850 < 26.5. The color selection incompleteness limits are roughly the same when applied 
to our model and the empirical Monte Carlo simulations preformed by Giavalisco et al. 
(2003b). We select 67% of all galaxies and 62% of stellar mass in that redshift range, down 
to ^850 < 26.5. This lends further credence to apphcability of the 5435-dropout color selection 
technique outlined in Giavalisco et al. (2003b), which was designed based on simpler, more 
empirical models of galaxy SEDs at this epoch, and has not, as yet, been extensively verified 
with spectroscopic data. 

We found a relatively good agreement between the model-derived and observed CDF-S 
i?435-dropout galaxy 2775 — colors. The ^775 — 2:850 colors showed much less agreement, 
however, we attribute the discrepancy to potentially inadequate dust recipe employed in the 
model. In general, we found that the stellar age, mass, and the star formation properties 
of the 2; ~ 4 sample were similar to the 2; ~ 3 sample studied by SPF2001. We found a 
model-derived ~ 40% stellar mass build-up for the UV rest-frame L* galaxies between the 
two epochs. In addition we found that if we look at all galaxies predicted by our model, 
spanning 3.44 < z < 4.12 redshift range, then at least 50% of the stellar mass contained in 
objects with Log[Mstar] > 9.26 (M© h~^) is missed due to the 2;85o < 26.5 observational limit. 
If our model is correct, and it does show reasonable agreement with respect to colors (modulo 
dust), number counts, and other physical quantities, this result has significant implications 
for the completeness of the optically-selected surveys at this epoch. 
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Space Telescope Science Institute, which is operated by the Association of Universities for 
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Support for this work, part of the Space Infrared Telescope Facility (SIRTF) Legacy 
Science Program, was provided by NASA through Contract Number 1224666 issued by the 
Jet Propulsion Laboratory, California Institute of Technology under NASA contract 1407. 



REFERENCES 

Baugh, C. M., Cole, S., Frenk, C. S., & Lacey, C. G. 1998, ApJ, 498, 504 



-8- 



Calzetti, D. 1997, AJ, 113, 162 

Chariot, S., Fall, S. M. 2000, ApJ, 539, 718C 

Devriendt, J. E. G., Guiderdoni, B., & Sadat, R. 1999, A&A, 350, 381D 

Pranx, M. et al. 2003, ApJ, 587, L79 

Giavalisco, M. et al. 2003a, this ApJL issue 

Giavalisco, M. et al. 2003b, this ApJL issue 

Giavalisco, M., Steidel, C., & Macchetto, D. 1996, ApJ, 470, 189 

Madau, P., Ferguson, H., Dickinson, M., Giavahsco, M., Steidel, C., & Pruchter, A. 1996, 
MNRAS, 283, 1388 

Meurer, G., Heckman, T., & Calzetti, D. 1999, ApJ, 521, 64 

Oke, J. B. 1974, ApJS, 27, 21 

Papovich, C. et al. 2003, this ApJL issue 

Papovich, C. 2002, Ph. D. thesis. The Johns Hopkins University 
Papovich, C., Dickinson, M., & Ferguson, H. C. 2001, ApJ, 559, 620 

Shapley, A. E., Steidel. C. C., Adelberger, K. L., Dickinson, M., Giavahsco, M., & Pettini, 
M. 2001, ApJ, 562, 95 

Somerville, R. S. 2002, ApJ, 572, L23 

Somerville, R. S., Primack, J. R., & Faber 2001, ApJ, 310, 1087 
Somerville, R. S., & Primack, J. R. 1999, MNRAS, 310, 1087 

Steidel, C., Adelberger, K. A., Giavahsco, M., Dickinson, M., & Pettini, M. 1999, ApJ, 519, 
1 

Steidel, C., Giavalisco, M., Dickinson, M., & Adelberger, K. 1996a, AJ, 112, 352 

Wechslcr, R. H., Somerville, S. S., Bullock, J. S., Kolatt, T. S., Primack, J. R., Blumenthal, 
G. R., & Dekel, A. 2001, ApJ, 554, 85 



This preprint was prepared with the AAS macros v5.0. 



-9- 




Fig. 1. — The ^775 — Kg versus ^775 — 2:850 colors of the z ~ 4 color-selected galaxies. Panel 
(a) shows the 136 observed S435-dropout galaxies filled- diamond from the CDF-S. Panel (b) 
shows a Hess diagram of the 'CoUisional-Star burst' model galaxies with artificial observa- 
tional scatter folded in, and selected with the same color criteria and magnitude limits that 
were applied to the observed sample. Contours have been superimposed to help guide the 
eye; they represent the 30*^, 50*^, 70*^, and 90*'' percentiles. In addition, histograms are 
included for both the model {gray solid-line) and observed {black dash-line) galaxy colors. 
This figure illustrates the relative agreement in the 2775 — Kg color distributions for the 
model-derived and observed 5435-dropout galaxies, while at the same time illustrating an 
appreciable mismatch in the ^775 — 2:850 colors. 




Fig. 2. — The stellar masses of the color-selected model galaxies. The filled- diamond symbols 
show stellar masses of the individual color-selected model galaxies vs. their corresponding 
^850 magnitudes. The histogram shows the projected distributions for the same color-selected 
model galaxies with an imposed 2:850 < 26.5 magnitude limit (shown as a gray dash-line). 
The predicted masses are two magnitudes lower than the stellar masses of the present day 
spirals and ellipticals. 
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Fig. 3. — This figure illustrates the stellar mass distribution of all mo del- derived galaxies 
spanning the redshift range of 3.44 < z < 4.12 and with Log[Mstar] > 9.26 (M0 h^^) 
(the median value from Figure 2). Galaxies were binned into 0.25 magnitude intervals and 
weighted by their corresponding stellar mass. A gray dash-line is included to delineate 
our observational magnitude limit. This figure shows the amount of model-predicted mass 
potentially missed by the current optical surveys. 



